Quantum state reconstruction with binary detectors 



D. Mogilevtse\|3 

Institute of Physics, Belarus National Academy of Sciences, F.Skarina Ave. 68, Minsk 220072 Belarus; 
Departamento de Fisica, Universidade Federal de Alagoas Cidade Universitdria, 57072-970, Maceio, AL, Brazil 

I present a simple and robust method of quantum state reconstruction using non-ideal detectors 
able to distinguish only between presence and absence of photons. Using the scheme, one is able 
to determine a value of Wigner function in any given point on the phase plane using expectation- 
maximization estimation technique. 
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Development of effective and robust methods of quan- 
tum state reconstruction is a task of crucial importance 
for quantum optics and informatics. One needs such 
methods to verify the preparation of states, to analyze 
changes occurring in the process of dynamics and to in- 
fer information about processes causing such a dynamics, 
to estimate an influence of decoherence and noise-induced 
errors, to improve measurement procedures and charac- 
terize quantum devices. 

For existing schemes of quantum state reconstruction 
losses are the major obstacle. In the real experiment they 
are unavoidable; detectors which one has to use to collect 
the set of data necessary for the reconstruction are not 
ideal. Presence of losses poses a limit on the possibility 
of reconstruction. For example, in quantum tomography 
P, 0] , which is up to date is a the most advanced and 
successfully realized reconstruction method, an efficiency 
of detection should exceed 50% to make possible an in- 
ference of the quantum state from the collected data. 
However, the very presence of losses can be turned to 
advantage and used for the reconstruction purposes. 

In 1998 in the work Q it was predicted, that non-ideal 
binary detectors can be used for complete reconstruction 
of a quantum state. A detector able to distinguish only 
between presence and absence of photons is able also to 
provide sufficient data for the reconstruction. This de- 
tector must be non-ideal, since the ideal binary detectors 
measures only the probability to find the signal in the 
vacuum state. 

To perform the reconstruction one needs a set of probe 
states mixed via a beam-splitter with a signal state. For 
the probe coherent states were suggested. Wheir the 
probe was assumed to be the vacuum, the procedure gives 
an information sufficient for inference of a photon num- 
ber distribution of the quantum state. In Inference of 
the photon number distribution was discussed in works 
P|. This scheme was implemented experimentally to re- 
alize a multichannel fiber loop detector fi . Very recently 
it was developed further by implementing the maximal 
likelihood estimation realized with help of expectation- 



maximization (EM) algorithm |y| , and demonstrated ex- 
perimentally !^ . The reconstruction procedure with 
help of EM algorithm was shown to be robust with re- 
spect to imperfections of the measurement procedure 
such as, for example, fluctuations in values of detector's 
efficiencies. In difference with the quantum tomography 
reconstruction scheme 0, , such a procedure does not 
impose lower limits on detector's efficiency and requires 
quite a modest number of measurements to achieve a 
good accuracy of the reconstruction. 

Here we demonstrate how to reconstruct a quantum 
state using sets of binary detectors with different efficien- 
cies. Let us consider a following simple set-up: the signal 
state (described by the density matrix p) is mixed on a 
beam-splitter with the probe coherent state Then 
the probability p to have no counts simultaneously on 
two detectors is measured (as it can be seen later, in 
fact, it is possible to use only one detector for the recon- 
struction) . 

According to Mandel's formula, this probability is 

p = {: exp{-Vcc) Vdd}d} (1) 

where Vc^ Vd are efficiencies of the first and second detec- 
tors; c\c and d\d are creation and annihilation opera- 
tors of output modes and :: denoted the normal ordering. 
For simplicity we assume here, that there is no 'dark cur- 
rent', and in absence of the signal detectors produce no 
clicks. Let us assume, that the beam-splitter transforms 
input modes a and b in the following way 

c = acos(Q;) -f 6sin(Q;), d = 6 cos (a) — a sin (a) . (2) 

Then averaging over the probe mode 6, from Eqs. 
and ||2Jl one obtains 

p = eyTr{: exp {-D{a) -f 7*)(a + 7)} : p}, (3) 

where 

V ^ Uc cos^(a) + Ud sin^(Q!), 
7 = I3{vd ~ Vc) cos(a) sa\{a)/u, (4) 
y = -\l3\'^VcVdSiv?{2a)/v. 

Finally, from Eq. ^ one obtains 
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p^eyY^{l-Dr{n\D\^)pD{^)\n), (5) 
n=0 
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where -D(7) = exp{7a''^ —7*0} is a coherent shift oper- 
ator, and \n) denotes a Fock state of the signal mode 
a. 

Essence of the reconstruction procedure is in measure- 
ment of p for different values of 9 and fixed value of the 
parameter 7. Let us, for example, assume detectors ef- 
ficiencies z/c, i^d to be constant, and I'c 7^ i^d (one of ef- 
ficiencies can be set to zero; only one detector might be 
used in the scheme). Then for arbitrary 7 let us mix 
the signal state with the probe coherent state having the 
amplitude 

^ (z^d - I'c) sin(2aj ) 

and measure a set of probabilities p for different values of 
the beam-splitter rotation angle aj . Then we have linear 
positive inverse problem of finding quantities 



i?„(7) = {n\D^j)pD{j)\n 



(7) 



which could be solved by the EM iterative algorithm sim- 
ilar to the one used in works 0, ll| for the reconstruction 
of diagonal elements of the signal state density matrix. 
Besides, Rn{0) are diagonal elements of the density ma- 
trix of the signal. 

The EM algorithm for the suggested scheme of the re- 
construction is as follows. We assume the signal density 
matrix to be finite in Fock-state basis being N x N, and 
the number of different values of aj is M > N. Then we 

chose an initial set of Rn\'y) > 0, Vn, and implement 
the following iterative procedure 

J=0 



(i-^,)"pr 



(8) 



where p^^^ is the experimentally measured frequency of 
having no clicks on both detectors for a given aj, and 

Pj'^'' is the left-hand side of Eq. (pj calculated using the 
result of fc-th iteration. The weights 

N-l 



The procedure (jSJ guarantees positiveness and unit sum 
of the reconstructed i?„(7). 

Finally, having reconstructed quantities i?„(7), it is 
straightforward to find a value of the Wigner function at 
the point 7 [Tol |: 



N-l 

m7*,7) = -E(-l)"^"(^)- 

TT ^ — ' 

n=0 



(9) 



The scheme can be even more simplified if we set the 
efficiency of the second detector to zero Vd = 0- Then 
p is simply the probability to have no clicks on a single 



detector. We have also y = 0, and the parameter 7 does 
not depend on the efficiency: 



7 = — /3 tan(a). 



(10) 



Thus, one can determine Rn{l) by varying the effi- 
ciency Vc a-nd keeping a constant without even changing 
the amplitude of the probe coherent state /?. 

In Figure n the reconstruction of the Wigner function 
of the signal coherent state with help of the one-detector 
version of the method is illustrated. 




FIG. 1: Reconstruction of the signal coherent state a with 
the amplitude a = 1. In Figure (a) the reconstructed Wigner 
function is shown; in figure (b) the difference between the 
exact Wigner function and the Wigner function of the trun- 
cated state is shown. Figure (c) shows difference between the 
exact and the reconstructed Wigner functions; in Figure (d) 
the variance (7(7,7*) is depicted. For all pictures Nr = 10* 
measurements were used for each point on the phase plane 
and Nit ~ 10'^ iteration of the EM algorithm. The state was 
truncated with A'^ = 12; 30 different values of the detector 
efficiency were taken; they were distributed homogeneously 
in the interval [0.1, 0.9]. 

Starting from the miiform distribution Rn\'^) = 1/N , 
very good results of the reconstruction was achieved with 
only a 10'^ iterations of the EM algorithm and 10^ mea- 
surements for each point on the phase plane. As it was 
mentioned in the work 7], the choice of Rn\'^) 7^ 
was indeed not infiuencing much the convergence for any 
given point. However, for different points on the phase 
plane the rate of convergence might differ strongly. In 
region of more rapid change of the Wigner function one 
needs more iterations and more measurements to achieve 
the same precision (as it can bee seen in Figure^d); here 
the variance is smaller near the peak of 1^(7,7*)). An 
explanation can be easily found from the formula (|5J| : in 
the region of rapid change one needs to find with high 
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precision several comparable i?„(7), whereas, for exam- 
ple, the behavior of the Wigner function near 7 = a is 
defined mostly by i?o(7)- Also, the precision is influenced 
by the truncation number TV. Increasing the region on 
the phase plane, where the Wigner function is to be esti- 
mated, one needs also to increase N. In Figure ^b) one 
can see the difference between the exact and truncated 
Wigner functions. As a consequence (as Figuresmc) and 
(d) show) , in the regions when the truncation error is sig- 
nificant errors of the reconstruction procedure are also 
increased. 

For illustration of how the total error of the reconstruc- 
tion propagates, we use the average distance between val- 
ues of the exact and the reconstructed Wigner functions 



'^^ = ^ E \We.act{l*,l) - 1^(7*, 7)1, 

V7 



(11) 



where the summation taken over all points on the phase 
plane were the estimation was made; 
number of points on the phase plane 



Np stands for the 
It can be seen in 



FIG. 2: The propagation of the reconstruction error ^VK lllll 
for different numbers of iterations Nit in dependence on the 
number of experiment runs Nr; for all curves Np = 2500 
and the following region on the phase plane was taken 
Re{'y),Im{j) G [—1.2,2.5]. Other parameters are as in Fig- 
ure 

Figure 12 that for smaller number of iterations Nn an in- 
crease of Nr leads to quicker convergence for small num- 
ber of measurements; after that the error AW decreases 
very slowly with increasing of Nr. An increase in the 
number of iteration leads to much slower convergence 
for small Nr. However, with increasing of Nr an accu- 
racy improve more rapidly; the error goes below values 
achieved for smaller number of iterations. Generally, Fig- 
ure |2] confirms an observation made in the work p|: for 
performing the reconstruction procedure it is reasonable 
to use a number of iterations close to the number of mea- 
surements, since for Nr ^ Nit an accuracy improves very 



slowly, and too large Nn might even lead to increasing of 
AW. 

From the practical point of view, to perform the re- 
construction it is reasonable first to estimate a photon 
number distribution (to find the set of i?„(0)). This will 
provide a clue for estimating the region of the plane suffi- 
cient for the reconstruction, and also the truncation num- 
ber necessary for the purpose. 



(a) 



(b) 
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FIG. 3: Reconstruction of the squeezed coherent state Hl.'ifl 
with tanh(r) — 0.5. In Figure (a) the reconstructed Wigner 
function is shown; in Figure (b) one can see Pmn obtained 
with help of Eq. 1121 . In Figure (c) the difference between 
exact and reconstructed pmn is shown; Figure (d) shows the 
variance amn. For all Figures the Wigner function was found 
in Np — 2500 points; the following region of the phase plane 
was used Re{'y) G [—1.1], Im{'y) G [—3,3]. Other parameters 
are as in Figure Q 

It is possible to infer elements pmn of the signal state 
density matrix in the Fock state basis using the recon- 
structed Wigner function in the following way [Tlj ]: 



d'7(-l)"W/(7*,7)^m«(27), 



(12) 



where 



A„„(27) = exp{-2|7|2}V';^ X 

min(m,n) , , , , , ™ i 

^ (27)"-'(-27*)™^' 
n(m- /)!(n- ?)! ' 

In Figure 13 an example of matrix elements pn 
squeezed vacuum state 

|r) =exp{-r2(at2-a2)/2} 



of the 



(13) 



is demonstrated. One can see that even for modest num- 
ber of points on the phase plane (50 points along each 
axis) and measurements (10^ per point) the accuracy of 
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the reconstruction of pmn is remarkable. The truncation 
error does not influence this elements much, because the 
function D,nn{2j) is small in the regions where the trun- 
cation error strongly influences the reconstructed Wigner 
function. One can mention here, that for reconstruction 
of Pmn from the quantities Rnil) one does not need to 
find the Wigner function in the whole region required for 
the integral J (P'jW{'j,j*) to be close to unity. It is suf- 
ficient to find Rnil) on N circles on the phase plane [T^ . 
However, this approach leads to non-positive problem of 
inferring pmn from the set of Rnij)- 

To conclude, we have suggested and discussed a sim- 



ple and robust method of the reconstruction of the quan- 
tum state of light. For the method one needs to use 
a binary detector, a coherent state for the probe and a 
beam-splitter. The reconstruction problem is linear and 
positive, and solved with help of the fast and efhcient EM 
algorithm of the maximal likelihood estimation. With 
help of the method a value of the Wigner function of the 
signal state can be found in any point on the phase plane. 
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